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Abstract 

We consider an emulsion whose droplets contain a trapped species (insoluble 
in the continuous phase), and study the emulsion's stability against coarsening 
via Lifshitz-Slyozov dynamics (Ostwald Ripening). Extending an earlier treat- 
ment by Kabalnov et &l(Colloids and Surfaces, 24 (1987), 19-32), we derive a 
general condition on the mean initial droplet volume which ensures stability, 
even when arbitrary polydispersity is present in both size and composition of 
the initial droplets. We distinguish "nucleated" coarsening, which requires either 
fluctuations about the mean field equations or a tail in the initial droplet size 
distribution, from "spinodal" coarsening in which a typical droplet is locally un- 
stable. A weaker condition for stability, previously suggested by Kabalnov et al., 
is sufficient only to prevent "spinodal" coarsening and is best viewed as a con- 
dition for met ast ability. The coarsening of unstable emulsions is considered, and 
shown at long times to resemble that of ordinary emulsions (with no trapped 
species), but with a reduced value of the initial volume fraction of dispersed 
phase. We discuss the physical principles relevant to the stability of emulsions 
with trapped species, describing how these may be exploited to restabilise par- 
tially coarsened emulsions and to "shrink" previously formed emulsion droplets 
to form "miniemulsions" . 

1 Introduction 

A phenomenological model for the growth of small droplets or precipitates from a su- 
persaturated phase, was first investigated by Lifshitz and Slyozov|I|] in 1961. They con- 
sidered a two-phase system consisting of a vanishingly small volume fraction of droplets 
of a dispersed phase, evenly distributed throughout a second, continuous phase. The 
model is extremely general and forms an explanation for a number of different physical 
phenomena. Among these is the coarsening of emulsion droplets made of a liquid of 
low but finite solubility in the surrounding continuous phase 0, [3], [|, [5], |6|. The Lifshitz- 
Slyozov model is based on a mean-field description of the evaporation/condensation 
process (Ostwald ripening), which ignores the contributions of coalescence (droplet 
fusion). Depending upon the nature of the local interactions between droplets, their 
mobility, and the solubility of the dispersed phase within the continuous one, this can 
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be a good approximate description of the coarsening of emulsion droplets. The effects 
of coalescence are ignored in this paper. The absence of coalescence may require the 
presence of surfactant to stabilize the emulsion droplets (for example by imparting a 
surface charge). We do not treat surfactant dynamics explicitly, assuming in effect that 
the surfactant transports rapidly through the continuous medium to ensure a constant 
surface tension of droplets at all times. 

Recent work|7J, || |9], has considered how the process differs when there are 
(disregarding surfactant) three species present [0, |3|, §]. In particular, Kabalnov et al.fT] 
considered a limiting situation where the third species is completely immiscible with 
the continuous phase, and is therefore effectively trapped within droplets. Kabalnov et 
al. gave quantitative criteria for the stabilisation of an emulsion by trapped species 
under the condition that the initial droplets were narrowly distributed in both size and 



composition. Related ideas to these can be found in the literature on meteorology fll 



12, [H^, 14, 15], where (for example) the stability of sea-mists and/or smog is in part 
attributed to the inhibition of coarsening by the entrapment of salt or other non volatile 
species. 

In this paper, we consider the effects of trapping a third species within emulsion 
droplets in greater detail. (Although we speak only of emulsions, several of our ideas 
relating to the stabilisation of droplets and precipitates to coarsening may be useful in 
other contexts |T5], [IB, 17 1.) We reexamine the condition for stability proposed in Ref.]7J, 



and conclude that this is, for monodisperse droplets, the requirement for metastability, 
rather than full stability. We give a more stringent condition for full stability, which 
allows for fluctuations about the mean-field growth rate of droplets of a given size. 
Our stability condition is then extended to the case of an arbitrary initial distribution 
of droplet sizes and compositions. We note that this stability condition generalises 
to nonideal equations of state, such as would be needed to describe trapped salt, for 
example. The classification of stable and metastable regimes is, for emulsions, somewhat 
different from the one which applies to a single droplet surrounded by vapour, as was 
studied (with a trapped species present) in Ref. M. 

These arguments, which are described in Sections 2-4, are both general and rigourous. 
Indeed, in the absence of coalescence processes, the problem of the stability of dilute 
emulsions with trapped species is found to be analogous to the finding of phase equi- 
librium in multicomponent mixtures. The only difference is that some extra thermody- 
namic variables are required, to describe the fixed number of droplets and their trapped 
contents. These variables would not be conserved if either coalescence were to occur, or 
the trapped species were able to diffuse through the continuous phase at a finite rate. 
(The latter case is briefly discussed, at the end of Section 4.) But if these processes are 
slow enough, our methods allow the state of the system, after the remaining degrees of 
freedom have attained equilibrium, to be determined. 

In Section 5, we consider the effect of a trapped phase on the kinetics of coarsening, 
in a system which does not obey the criterion for stability derived earlier. Approximate 
solutions for the droplet size distributions are given, and the physically motivated 
approximations upon which they are based are explained. Our analysis supports the 
view of Ref. 0, that in the unstable (coarsening) growth regime, the long time dynamics 
is essentially the same as given by the Lifshitz Slyozov description. 

In Section 6, several implications and possible applications of these results are 
considered, including possibly novel methods for forming stable emulsions of small 
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droplet size, and ways of arresting or reversing coarsening after it has already begun. 
Section 7 gives a brief conclusion and outlook for further work. 



2 Kinetics of Droplet Growth 

We consider a vanishingly small volume fraction of droplets and assume that droplets 
are in equilibrium with their local spatial environments. We also assume that the en- 
semble of droplets interacts via diffusion of the dispersed-phase species through the 
continuous medium, and that this diffusion is rate-limiting (it is much slower than the 
rate at which the dispersed-phase species may be absorbed into a droplet from imme- 
diately outside). Each droplet's growth rate is then determined by the difference in 
concentration of dispersed-phase species immediately outside its surface and the (aver- 
age) far-field concentration within the continuous phase. These assumptions reduce the 
calculation to (i) a single body equilibrium problem, to find the mean growth rate of a 
droplet in terms of its size and the ambient supersaturation; and then (ii) a many-body 
problem to determine the droplet size distribution as a function of time. Since in stage 
(ii) each droplet interacts with the others only through the ambient supersaturation, 
the problem is tractable. 

The method and assumptions just described are basically the same as those first 
used by Lifshitz and Slyozov||T| to obtain the average rate of droplet growth and an 
asymptotic solution for the droplet size distribution. Although other methods have 
been used to correct and improve those used by Lifshitz and Slyozov, the resultant 
corrections are small and the resulting asymptotic solutions essentially the same. For 
a clear exposition of the main ideas see Bray[Q. A review of recent theoretical and 



experimental work is given by Voorhees[|19|j, and includes an overview of theoretical 
work on systems at finite droplet volume fractions. 

2.1 Local Equilibrium of a Droplet 

We treat all components of the emulsion as incompressible fluids, in which case the 
Helmholtz free energy of a dispersed-phase drop is given by 

F d = F b + aATtR 2 (1) 

where Ff, is the free energy of the same amount of bulk liquid, a is the surface tension, 
and R is the droplets radius. This results in a chemical potential given by 

2av b 

Pd = Pb + (2) 

where v b is the volume of a single molecule, and \x b the chemical potential, in a bulk 
liquid of the dispersed phase. 

Some of the dispersed-phase species does not reside in droplets but is instead molec- 
ularly dissolved in the continuous phase at concentration C. The chemical potential of 
such molecules is 

/i c = /i c o + k B T In C (3) 

where /i c o is a reference value. (This assumes that the molecular solution is dilute 
enough to behave ideally.) At equilibrium = \i c so, immediately above a droplet 
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surface, the concentration C(R) obeys 



k B THC(R)) = ^ + »b-ti* (4) 

At a flat interface with radius of curvature R — > oo, we have k B T ln(C (oo)) = fib — £t c o, 
so that we may write 

C(R) = C(oo) exp (^L) (5) 
where /x c o has been eliminated and 

A 2(TVb fd\ 

Afi = fi d - fi b = — - (6) 
ri 

For small values of Afi/k B T, as are generally maintained during the late stages of the 
coarsening process, we may expand Eq.||] to obtain 

In fact, for typical emulsion systems near room temperature with interfacial tensions 
< 10^ 1 Nm^ 1 and initial droplets ~ 0.1/im or larger in size (as would arise for emulsions 
made by mechanical agitation) the approximation Afx/k B T <C 1 is likely to be valid 
throughout the coarsening process. 



2.2 Rate of Droplet Growth 

Under the conditions studied here c, the average concentration of dispersed-phase 
species dissolved in the continuous phase, will itself be close to C(oo). Hence the re- 
duced supersaturation of dispersed-phase species, which we define as 

e = (c-C(oo))/C(oo) (8) 

will be small. This means that droplet sizes change slowly on the time scale of relax- 
ation of the diffusion field: a steady-state approximation may then be made for the 
concentration profile around a droplet. This entails replacing the diffusion equation 
dc/dt = DV 2 c with V 2 c = 0. Imposing as boundary conditions c(R) = C(R) (the 
equilibrium value at the droplet surface) and c(oo) = c (the mean value far away), we 
obtain the steady-state profile 

c{r)=R (mz±y- c (9) 

for the concentration field c(r) at a distance r from the centre of a droplet of radius R. 

Since we assume that the process is diffusion-limited, the rate of droplet growth at 
the surface is given by the incident flux as 



dR_ vE) dc 
dt dr 



R 




(10) 
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Using the linearised form ([/]) for C(R) we obtain 



dR _ Dv b C(oo) / 2av b \ 

tZt ~ i2 V k B TR) [ } 

If we define reduced variables R' = Rk B T ' j '(2av b ), f = tDC(oo)k 2 B T 2 j '4v b a 2 and and a 
reduced growth rate U(R', e) = dR'/dt' then Eq.|ll] simplifies to|T8| 

This shows the balance between condensation (first term) and the evaporation driven 
by the Laplace pressure effect (second term). Note that no finite droplet size can be 
stable against evaporation at long times, when the supersaturation e tends to zero. 
Accordingly, emulsion droplets of a single species with finite diffusivity will always 
coarsen. 



3 The Effect of a Trapped Species 

Molecules of a third species which are entirely immiscible with the continuous phase, 
are now considered to be present within droplets. Such molecules are effectively trapped. 
This has the immediate consequence that no droplet of the dispersed phase can ever 
entirely evaporate. So the trapped phase results in the additional thermodynamic con- 
straint that the total number of droplets remains constant in time, with a value deter- 
mined by the initial conditions. 

We assume for simplicity that the trapped species may be treated as dilute within 
each droplet. (This assumption is relaxed in Section ^L7\ below). In this case an addi- 
tional term which corresponds physically to the osmotic pressure of the trapped phase, 
will contribute to the droplets chemical potential. Indeed it is easily confirmed]?] that 
Eq|] is modified to 

' 2a r}k B T 



A/i 



v b (13) 



R (4/3tt)R 3 _ 

where rj is the number of trapped particles in the droplet. Thus, the osmotic pressure 
of the trapped species (second term) competes directly with the Laplace pressure (first 
term); the latter favours fewer, larger droplets whereas the osmotic pressure favours a 
uniform droplet size (at least if all droplets have the same rj). 

Note that the second term diverges as R — > 0; although the trapped species cannot 
be treated as dilute in this limit, the above formula is already enough to prevent a 
droplet from ever evaporating completely. Using Eq.|l^ to determine a growth rate 
dR/dt we find that Eq.ll becomes 

dR _ Dv b C(oc) ( _ 2av b r)v b \ 

dt~ R \ k B TR (47r/3)i? 3 y 1 ' 

When rj = 0, a characteristic length scale R t can be defined by 

ek B T 
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which gives a critical size above which droplets will grow in size and below which 
droplets will dissolve. When 77 > we may define another characteristic length scale 
Rb by that at which the osmotic and Laplace pressures are in balance, which gives 

Hence the growth rate may be written in the more informative manner 

dR Dv 2 b C(oo)2a 



dt Rk B T 




(17) 



As shown in figure [I], when rj > a second zero at R ~ Rb is introduced into the growth 
rate in addition to the original one. The stability of the new fixed point at R ~ Rb 
suggests the possibility of forming stable emulsions; but since the curves depend on e, 
which is time-dependent and need not even tend to zero at long times, the ultimate 
behaviour of the system is unclear. This is considered further in the following Section. 
It is possible to define a growth velocity U(R', e) in terms of dimensionless variables, 

as 

= £-5> + iF (18) 

where R',t' are as defined in Section |2| above, and i?'J = 3r]k B T 3 Vb/32TT<j 3 . This form 
of the growth rate is used in Section |5| when considering the dynamics of unstable 
distributions which coarsen. 



4 Formation of Stable Emulsions 

In discussing stability criteria, we assume a monodisperse initial condition with an 
ideal equation of state for the trapped species within the droplets. These assumptions 
will be relaxed in Sections £L5| and [L7| below. 

4.1 Equilibrium at Fixed Droplet Number 

Consider an initially monodisperse emulsion in which the number of trapped particles 
77 is identical in each droplet. Since the emulsion droplets are treated as macroscopic 
objects their entropy of mixing is negligible, and may be ignored. The total free energy 
density of the system, at some later time, may be then written 

F = J dVn(V)f(V,r]) (19) 

where n(V) is the number density of droplets of volume V and f(V, 77) is the free energy 
of such a droplet containing 77 trapped molecules. For an ideal trapped species, f(V,r)) 
takes the form 

f(V, rj) = n b V/v b + 47ra(3y/47r) 2/3 - r]k B T In V (20) 

where a constant has been suppressed. The curves with and without trapped species 
are shown in figure ||. 
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Figure 1: A schematic diagram representing how the presence of a trapped phase in- 
troduces a stable droplet size. Droplets of a given size move in the direction shown 
by the arrows. Without trapped species there is one unstable fixed point at Rs- With 
a trapped species a new, stable, fixed point appears at Rb due to the competition 
between Laplace and osmotic pressure. If there are trapped species and e is negative 
then there is again only one fixed point, but it is the stable one at Rb- Formulae for 
Rs, Rb appear in Section 4.2 below. 

The free energy density F is to be minimised subject to the constraints that 
/ dV n(V) = n (the initial number of droplets is held constant), and / dVVn(V) = 
4> = n V, representing a fixed total volume fraction <ft of dispersed phase |2(J. The lat- 
ter ignores a contribution from the dispersed-phase species that is solubilized in the 
continuous phase; for emulsions of low solubility, prepared by mechanical dispersion 
(rather than by quenching from a homogeneous mixture at high temperatures) the 
latter is always negligible. Note that the linear term, HbV/vb in /(V) is irrelevant to 
this procedure and, for clarity, we have subtracted it from the curves shown in Fig.|] 
and similar curves appearing below. 

These rules are precisely analogous to those for constructing the equilibrium state 
of a system whose free energy /(A) depends on a composition variable A. Such a system 
can separate into volumes v(X) of phases with different compositions A; in that case, 
F = Z)a' u (A)/(A) is minimised subject to the constraints J2x v W — v o (the total 
volume of the system is fixed) and J2x v W^ — v o^o (the total amount of species A is 
conserved). Here the subscripts describe a hypothetical homogeneous state. 

In this analogy, each value of the droplet size V corresponds to a "phase" of the 
system, and n(V) corresponds to the "volume" of such a phase present in the final 
state. A monodisperse emulsion represents a "single phase". The minimisation of F 
is therefore exactly as one would perform to find phase equilibrium in a binary fluid 
with the function /(A) replacing f(V,rj): the usual construction for this is to seek 
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Figure 2: Comparison of curves of f (V) — / v h without (left) and with (right) trapped 
species. 



common tangencies whereby F can be lowered by phase separation. (In this analogy, 
phase separation corresponds to the formation of droplets of more than one size.) The 
volumes v(X) are determined by the lever rale^H]], and, in principle, the same rule 



would apply here to calculating n(V) under conditions where more than one droplet 
size was present in equilibrium. Notice that the correspondence works only for the 
curve f(V, 77); no similar construction applies to f(R, 77), which is the form more usually 
considered in the literature J7|, ||. Also, note that any prediction of a "single phase" 
(i.e., a monodisperse emulsion) is, in principle, subject to a small spreading of the size 
distribution arising from the entropy of mixing of the emulsion droplets themselves 
(rather than of the contents within one droplet), which we have neglected. 



4.2 Stability Criterion for Monodisperse Emulsions 

According to the above argument, if 77 is the same for all droplets (as will usually be 
nearly true if the initial droplet distribution is monodisperse) the equilibrium state of 
the system at fixed droplet number can be found by inspection of the f(V, rf) curve. A 
monodisperse emulsion can be stable only if the corresponding V lies in a part of the 
curve of positive curvature; if this is not the case then the free energy may be reduced 
(at fixed total number of droplets and fixed 77 in each drop) by the monodisperse distri- 
bution becoming polydisperse, and coarsening will occur. However, though necessary 
for stability, positive curvature is not sufficient to prevent coarsening, since even when 
the curvature is positive at a point V = Vq, it may be possible to find a lower free 
energy by constructing a common tangent on the f(V) curve which lies below f(Vo). 
The shape of the f(V) curve shown in figure Q dictates that any such tangency must 
connect the point at V — > 00 to the absolute minimum of / — fibV/vf,; the latter arises 
at V = Vb(t}) and is discussed further below. 

The weaker of the two stability conditions (positive curvature) is precisely that 
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given by Kabalnov et alj7j: 



Hence Kabalnov et alp] reasoned that a sufficiently monodisperse initial distribution 
with droplet size Vq < Vg = 47ri?|/3 would be stable. However, the above argument 
shows this is actually a criterion for metastability. In other words, the criterion of 
Ref. actually separates initial distributions which coarsen immediately, from those 
where coarsening requires some fluctuation to bring it about. Since every droplet has 
a slightly different environment, such fluctuations are invariably present, though they 
are not included within the theory of Lifshitz and Slyozov. The kinetic mechanism for 
fluctuation-induced coarsening is discussed in the following Section [4.3| . The extent of 
the metastable region, and the likelihood of a suitable fluctuation actually occuring, 
are discussed further in Section [4.6| . 

In any case, the above thermodynamic analogy shows that full stability in fact arises 
for initially monodisperse emulsions if and only if 

This corresponds to the requirement that the initial state lies to the left of the absolute 
minimum in the function f(V,rj) — fibV/vb- Everywhere to the right of this minimum, 
a lower global free energy can be constructed by a common tangency between Vb and 
(formally) V = oo. This corresponds to a final, coarsened state in which a monodisperse 
emulsion of droplet size Vb = 4ttR b /3 and droplet number density no coexists with an 
"infinite droplet", which can be interpreted as a bulk volume of the dispersed-phase 
species. 

The subscripts B, S for Vb,s an d Rb,s can be taken to denote "balance" (between 
osmotic and Laplace pressures in a drop) and "stability" (in the sense of Kabalnov 
et al||). However, in view of the thermodynamic discussion, it might be better to 
interpret them as "binodal" and "spinodal". Indeed, depending upon the value of Vb, 
any monodisperse distribution lies in one of three possible regimes: see figure []. In 
regime I, the emulsion is fully stable under coarsening dynamics (though not, of course, 
under coalescence). In regime II it is metastable. In regime III it is locally unstable 
and will coarsen immediately. The metastable region II is both mathematically and 
physically analogous to that between the binodal and the spinodal lines governing 
phase coexistence in a binary fluid. 

4.3 Kinetic Interpretation 

The above thermodynamic argument is quite formal, so it is useful to interpret it in 
kinetic terms. Consider a monodisperse distribution which is in equilibrium with its 
"vapour" (i.e., the dissolved fraction of the dispersed-phase species). Let the radius of 
all but one droplet be (say) R < R$, but suppose a single larger droplet is present, 
of radius R. The growth of this droplet is determined by whether C(R,rj) is larger 
or smaller than c. Since the remaining droplets are in equilibrium, c will be equal to 
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Figure 3: The three different stability regimes. 



C(R ,rj). Now consider the concentration of dispersed-phase species C(R,r)) at the 
surface of the anomalously large droplet (see figure |j). 

Distributions in regime II have Rb < Ro < Rs and an average concentration of 
disperse phase c(Rq,t]) > C(oo,0). Hence if a perturbed droplet is sufficiently large, 
its surface concentration will be below the ambient level (see figure |j) and the droplet 
will grow as a result of diffusive flux through the continuous phase, at the expense of 
smaller droplets. Distributions in regime I (Rq < Rb) have c < C(oo,0) and there is 
no size of droplet larger than Rq for which sustained growth is possible. Such a droplet 
will instead redissolve to rejoin the equilibrium droplets at size R . 

Returning now to distributions in regime II (Rb < Ro < Rs), the critical size Rc 
above which a droplet will become unstable and start to grow is given by the larger 
root of dR/dt = 0, where dR/dt is given by Eq.[l7]. This gives 



R C = Re\l 



Rb 

~rZ 



R 



B 



R e 



(23) 



where R e is as in Eq.[L5] and depends upon the supersaturation e. Since the initial dis- 
tribution of droplets is considered to be monodisperse with size R and in equilibrium, 
then dR/dt\ft = and R e may be determined as 



Rr 



Rl 



Rn 



(24) 



This suggests that, for a monodisperse initial state, the nucleation time for the coars- 
ening process to begin (requiring nucleation of a droplet of size Rc) can become very 
large as Rq approaches Rb- On the other hand, nucleation can occur immediately if 
there is any slight tail to the initial size distribution, extending beyond Rc- 
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Figure 4: Variation in concentration of disperse phase at a droplets surface as a function 
of droplet volume. 

The distinction between nucleation due to the existence of a tail in the initial size 
distribution which extends to abnormally large droplets, and nucleation resulting from 
fluctuations in growth rates, is similar to that between heterogeneous and homogeneous 
nucleation in conventional phase equilibrium. In this analogy the tail of the droplet size 
distribution (which may be negligible for thermodynamic purposes), provides "nucle- 
ation centres" which allow coarsening to begin. 

4.4 Laplace and Osmotic Pressure Balance 

The absolute stability requirement, Rq < Rb (Eq.|22|), is equivalent to requiring that 
C(Ro,rj) < C(oo,0). However, we have for an ideal trapped species 

C(R,r]) = C (oo,0) (l + 

Identifying 2a /R with the Laplace pressure Hl and r/ksT / \Ati /3)R 3 with an osmotic 
pressure H OS m of the trapped species, we see that the absolute stability condition may 
be written simply as 

n L < n osm (26) 

or alternatively 2a /ksT < c t R where c t = rj/V is the concentration of trapped species 
in the dispersed phase. This means that if an emulsion is made from a dispersed phase 
of fixed Ct, a monodisperse emulsion will be fully stable only if the initial droplet size 
is sufficiently large. (The same applies to the metastability condition of Ref.[Q.) This 
may be a rather unintuitive result; indeed, at fixed number of trapped particles 77, as 
was used to discuss the f{V,r]) curves, stable emulsions arise only for small droplet 



2a 
~R 



yk B T 
{An/3)R 3 



(25) 
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sizes (Rq < Rb)- However, for a fixed composition q, the value of rj depends on i? ; 
and the important requirement is that the osmotic pressure (which inhibits coarsening) 
of an initial droplet exceeds its Laplace pressure (which drives it). At given ct this is 
true only for large enough droplets. 

4.5 Stability of Polydisperse Emulsions 

The criterion discussed above for full stability can be extended to the polydisperse case. 
Let us first consider the case where the initial state contains droplets with variable sizes 
Vq, but exactly the same 77. This system will again find its equilibrium state at fixed 
total number of droplets no and fixed volume fraction <p, which now obeys 

<p = n V (27) 

with V the mean initial droplet volume. By again ignoring the entropy of mixing of 
the droplets (treating them as macroscopic objects), we see that the thermodynamic 
arguments developed in Section [4.1| for equilibrium at fixed droplet number apply 
without modification to this type of polydispersity. Therefore, the ultimate behaviour 
is found simply by substituting V for V in our previous discussion. Thus the criterion 
for full stability in this case is V < Vb- It is also clear that in the stable regime the 
droplet distribution will evolve under the evaporation/condensation dynamics into a 
monodisperse one, whatever its initial polydispersity. This contrasts with the other 
regimes, where the final state will again have density no of monodisperse droplets, but 
of size Rb and in coexistence with an "infinite droplet" containing the excess amount 
of dispersed phase (n (V — Vb)) not residing in finite droplets. 

The situation when 77 is not the same for all droplets is more complex. To find 
a stability condition for this case, we consider first a more formal argument which 
reproduces the above result for a single 77. If the emulsion is unstable, we expect that 
at long times the size distribution will split into a "coarsening part" and a "stable 
part". As equilibrium is approached, the coarsening part has some average size V c 
which tends to infinity at long times. This requires that e — - > 0, which in turn means 
that the stable part of the distribution is necessarily monodisperse, and of droplet size 
V — > Vb- (This corresponds to the coexistence condition between finite and an infinite 
droplet mentioned previously.) 

Let no, n s and n c be respectively the number densities of all the droplets in the 
system, of those in the "stable part" and of those in the "coarsening part" of the distri- 
bution; clearly n = n s + n c . Then as the system tends to equilibrium, the conservation 
of (which we define to include the trapped species) requires 

n V = n s V B + n c V c (28) 

and hence 

n (V-V B )=n c (V c -V B ) (29) 

So if Vb > V, then > n c (V c — Vb); but since V c > Vb, then n c = and a coarsening 
part of the distribution cannot exist in this case. 

Now we consider the more general situation where there is polydispersity, not only 
in the initial droplet size, but also in the quantity rj of trapped species present in the 
initial droplets. Let n(rj), n s (r]) and n c {rj) be the number densities of droplets which 
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contain n trapped particles in the full distribution, and in its stable and coarsening 
parts respectively. As the system tends to equilibrium the conservation of <fi and no = 
J n(n) dn now requires 



n(n)Vdn = J n s (n)V B (n)dn + / n c {rj)V c {rj)dri (30) 

n s {rj) + n c (n) = n{rj) = n p(r}) (31) 

where p(n) is the probability of a given droplet having n trapped particles (which is 
time-independent). Hence 

n V = n J p(rj)VB{n)dr] . , 

+ Jn c (vWc(v)-V B {v))dr) ^ 

which may be rewritten as 

n (V - (V B ( V )) V ) = J n c (n)(V c (n) - V B {n))dn (33) 
SoifF< {V B {rj)) v , then 

0> / n c ( V )(V c ( V )-V B ( V ))drj (34) 



but since V c (r]) > V B (rj), then nc(n) = for all n, and no coarsening part of the 
distribution can exist. 

In summary, a condition which is sufficient to ensure the full stability of emulsions 
with an arbitrary initial distribution of sizes V and trapped species n is: 

V<(V B (rj)) v (35) 

where V = 4>/n and V B (rj) is as defined in Eq.||: V B (rj) = (4tt/3) 1/2 (rjk B T /2af /2 . 
The assumptions behind this result are: (i) the trapped species has zero solubility in 
the continuous phase and forms an ideal solution in each droplet; (ii) coalescence is 
strictly absent; and (iii) the solubility of the main dispersed-phase species is nonzero (so 
that diffusion can occur) but small enough that it contributes negligibly to 0. Subject 
to these, it is a rigorous result. It is easily established, using arguments that parallel 
those of Section |4.1| , that condition Eq.|35| is not only sufficient for full stability but also 



necessary, in the sense that any distribution which does not satisfy Eq.[35] can lower 
its free energy at fixed droplet number by a nucleated (if not a spinodal) coarsening 
process. As noted in Section [12] for the monodisperse case, for an emulsion which is 
only slightly unstable and whose size distribution does not have a tail extending to 
large droplets, the nucleation time before coarsening begins may be very long. If there 
is a tail, a long induction time is not expected. 

Notice that our rigorous condition Eq.|35| involves calculating (V B (r])) TI by averaging 
over the probability distribution of the trapped species. This quantity is not the same 
as V B ((n)); indeed for a > 1 one has the general inequality 

/ roo \ ct 

■r] a p(r])d7] > (J np(r])dr}j (36) 

Since V B ~ if I 2 (a = 3/2) we find that (V B (r])) v > V B ((rj}). Hence the approximate 
stability requirement V < V B ((rj)), based on the mean trapped particle number, un- 
derestimates the maximum initial droplet size. Accordingly this condition is sufficient, 
but not necessary, to ensure full stability. 
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4.6 The Metastable Regime (II) 

For a polydisperse system the condition for full stability is clear (see above), whereas 
that for metastability is less obvious. Even when all droplets have the same rj, metasta- 
bility may depend on the details of the initial droplet size distribution. Certainly, if this 
has no upper limit (i.e., there is a finite density of droplets above any given size), then 
the largest droplets will serve as nuclei for coarsening. As a qualitative rule, one can 
apply the argument of Section fO in a mean-field approximation, whereby an anoma- 



lously large droplet is considered to be exchanging material with a set of others, which 
for simplicity we treat as having a single size V. This enables a critical radius Rc to 
be estimated from Eq.[23| by replacing the length scale R € with 

As before, if the initial size distribution contains any droplets larger than R c , coars- 
ening can be expected. 

If 7] varies between droplets, things are still more complicated, since nucleation 
is likely to involve droplets of larger than average rj, as well as larger than average 
size. Although we still expect three regimes (fully stable, metastable and unstable) 
corresponding to those discussed in Section for the monodisperse case, the boundary 
between the metastable and unstable regimes may have a complicated dependence on 
the initial distribution of droplets and trapped species. This contrasts with the very 
simple criterion for full stability, Eq.[35] which applies for arbitrary initial conditions. 

Note that the dynamics we propose for nucleation in the metastable region is pecu- 
liar to emulsions. For example, Reiss and Koper|| discussed the equilibrium of a single 
drop within a uniform environment at fixed supersaturation. They concluded that for 
their problem nucleation dynamics were likely to be extremely slow, and deduced a 
criterion for stability corresponding to that given by Kabalnov et al. [0]. This does 
not contradict our own conclusions, because Reiss and Koper consider a single drop as 
opposed to a population of droplets. For an emulsion in regime II, only the largest drop 
present need exceed the nucleation threshold to initiate coarsening, and in many cases 
a large enough droplet would be present in the initial droplet size distribution. What 
is more, in a real emulsion with finite volume fraction, variations in the local environ- 
ments experienced by individual droplets will be a source of dynamical fluctuations 
which may far exceed the purely thermal fluctuations considered in Ref. ||. 

In summary, the stability conditions for a single droplet, and a population of 
droplets, are very different. We suspect that, even in a nominally monodisperse emul- 
sion, nucleation rates are generally not negligible, and hence that the only reliable 



criterion for stability is the one we have given in section 4.5 



4.7 Other Equations of State 

The thermodynamic arguments of Sections PO and fO generalise readily to an arbitrary 



equation of state for the trapped species. Indeed, for the case where rj is the same in 
all droplets, one need only replace Eq.|2D] with, 

f(V, 77) = fi b V/v b + 47Rx(3V/47r) 2 / 3 + f t (V, rj) (38) 
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where ft{V,rj) is the free energy of rj trapped particles in volume V and could include 
arbitrary interactions between these. If all droplets have the same rj, the condition 
for stability remains that V < Vb where Vsirj) is the droplet size corresponding to 
the absolute minimum of f(V,ri) — fi b V/vb- For such droplets, the osmotic pressure 
(—dft/dV)^ is again in exact balance with the Laplace pressure. The metastability 
criterion is again that f(V,rj) has positive curvature. For the case where rj is not the 
same for all droplets, the full stability criterion is again Eq.|35|. 

A complication arises if the interactions between the trapped species are attractive. 
Reiss and Koper|| pointed out that in this case a third, unstable fixed point can exist 
at a size smaller than Vb- They also correctly point out that such an attraction is 
likely to cause phase separation within droplets. The latter occurs whenever ft(V,rj) 
has negative curvature (with respect to V at fixed 77); the form of ft(V,rj) in Eq. 



38| must then be modified to reflect the internal phase separation, and any negative 
curvature regions will then be replaced by zero-curvature ones (corresponding to tie 
lines). In this case, because of the surface tension contribution in Eq.^, there can arise 
an additional minimum in f(V,rj) — LibV/vb which could lead, for example, to stable 
bidisperse emulsions for some range of initial conditions. The same can arise without 
intra-droplet phase separation, if ft(V, rj) has a small enough positive curvature for this 
to be outweighed by the negative contribution from the surface tension term in Eq.[38|. 
We leave a detailed discussion of these cases for future work. 

These results are sufficient, for example, to deal with the case of a distribution 
containing trapped salt (where ft(V, rj) can be approximated by, say, the Debye Hueckel 
equation of state p2|) which is of interest in meteorological as well as emulsion stability 



contexts [0, 0, |T|, [T|, [15 

In the case of repulsive interactions, at least, our stability condition Eq]35| is also 
easily generalised to a situation in which there is more than one trapped species. Letting 
rji be the number of particles of the ith species trapped within a droplet, then the 
condition for stability generalises to 

(39) 

where Vb(t]i, tj?, ■ ■ ■) is the droplet size corresponding to the absolute minimum of 
f(V,r)i,7)2, • • •) — fibV/vb- A proof of this generalized stability condition follows that 
given in Section [4.5| , with the single variable rj replaced by the list of variables r/x, 772, • • •■ 
So far, we have not considered explicitly the role of surfactant (which is usually 
present to prevent droplet coalescence in emulsions), tacitly assuming that this merely 
alters the constant value of a, the surface tension. For a surfactant that is insoluble in 
the continuous phase, this surface tension will itself be a function of droplet size. We 
do not treat this case further, but note that it could be included in Eq.[38] by replacing 
the term in a with a suitable "surface equation of state" for the trapped surfactant. 
The same applies to bending energy terms which could be significant for extremely 
small droplets. 

4.8 Behaviour when the "trapped" phase is partially soluble 

Throughout the above we have treated the trapped species as entirely confined in the 
emulsion droplets. In this Section we briefly consider what happens if this third species 
is very slightly soluble. (For simplicity we treat rj as the same for all droplets.) Our 
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previous classification into regimes I, II and III, though no longer strictly applicable, 
remains a guide to the resulting behaviour. The discussion that follows is related to 
that of Kabalnov et al.0. 

Recall that for entirely trapped species, the unstable regime (III) is characterised 
by the evolution of a bimodal distribution of droplet sizes, in which the larger droplets 
coarsen by the Lifshitz-Slyozov mechanism, while the smaller droplets adopt a size in 
equilibrium with the larger drops which approaches Vb as coarsening proceeds. (This 
scenario is examined in more detail in Section |5] below.) If the trapped phase is now 
made slightly soluble, the larger drops will, as before, coarsen at a rate determined 
by the transport of the majority (more soluble) dispersed phase. However, the small 
droplets that remain cannot now approach a limiting size, but will themselves evaporate 
at a much slower rate governed by the transport of the "trapped" species. Therefore, in 
the unstable regime, a two-stage coarsening is expected. In the fully stable regime (I), 
on the other hand, the first of these processes (rapid coarsening of the larger droplets) 
is switched off. The droplet size distribution has a single peak, which in the fully 
insoluble case will approach a delta function at the initial mean size. However with 
slight solubility coarsening will occur via the Lifshitz-Slyozov mechanism, but only at 
a slower rate controlled by the transport of the less soluble species. The behaviour in 
the metastable regime (II) is complex0, possibly characterised by a "crossover" from 
coarsening controlled by the less soluble component to coarsening controlled by the 
more soluble component, and we do not pursue it here. 

For a slightly soluble trapped species, the best prospect for stability is always to 
avoid the rapid coarsening process associated with the transport of the more soluble 
dispersed phase component. Since there is always some part of the size distribution 
that coarsens at the slower rate set by the trapped species, the requirement is satisfied 
so long as the size distribution remains single-peaked at all times. In Appendix |A| we 
show that a sufficient condition for such behaviour is: 

F(0) < V B (rj(0)) (40) 

where rj(t) is the average number of trapped species in a droplet at time t. Since 
Vb(j){Q)) < ^b(^(0)), this condition is less easily satisfied than our absolute stability 
condition, Eq.[35|, which applies for the case of entirely trapped species. Note, however, 
that although this condition is sufficient, we have been unable to find a necessary 
condition. In other words, it is possible that some systems will not satisfy Eq.^D] but 
will nonetheless coarsen only slowly. The derivation [^0| of Eq.[|0] requires Vb(t]) to 



increase faster than linearly with rj; this is valid for ideal mixtures and most other 
systems (perhaps excluding any in which interactions between the trapped species are 
attractive; see Section (4.7|). 



5 Coarsening Dynamics When a Trapped Species 
is Present 

5.1 The Physics of Coarsening 



We now consider conditions in which the requirement for stability, Eq.22, is not satis- 



fied. For simplicity we take the case where all droplets have the same number of trapped 
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particles r\. We also assume that any nucleation event required has taken place, and 
study the resulting coarsening process at late times. This entails small supersaturation 
e. Throughout this Section we work in the scaled variables R f , R' B and t' defined in 
sections [2] and |3|, although the primes are omitted for convenience. 

In the unstable mentioned in Section [| the equation for static equilibrium 

(zero growth) of a given droplet, U(R, e) = 0, has both a stable fixed point close to the 
balanced droplet size (R ~ Rb) and an unstable one at R ~ 1/e. The latter will lead 
to coarsening. The following arguments suggest that the presence of the additional, 
stable fixed point influences the coarsening dynamics only in a rather simple way. At 
long times, the main role of a population of stable droplets (R ~ Rb) is to effectively 
exclude a finite proportion of the dispersed phase from the coarsening process. The 
remaining part, whose volume fraction is effectively reduced, then coarsens almost 
normally. 

To see this, we first expand U(R, e) for small R — Rb and small e to obtain (in 
reduced units) 

U(R) = e/R B -(R- Rb) (^^) ( 41 ) 

with a root U — at 

R T (t) = R B + e(t)(^pj + 0(e 2 ) (42) 

At late times e will be small and therefore, according to Eq.^T], the size of any droplet 
in the neighbourhood of Rb relaxes exponentially toward Rr{t) with a fixed decay rate 
2/ R? B + 0(e). At long times, this process will be rapid compared to any coarsening of 
droplets of size R 3> Rb', hence any non-coarsening part of the droplet distribution 
comprises an effectively monodisperse population of radius Rr(t) obeying Eq.|42|. Put 
differently, Rr{t) is determined by e, because the non-coarsening droplets are in local 
equilibrium with the ambient supersaturation at all times. 

Now consider drops with size R ~ 1/e, in the neighbourhood of the second zero of 
U(R, e) (this is the unstable fixed point of the growth equation at given e). Expanding 
U(R, e) in R about R = 1/e we find 

U(R) = (R - l/e)e 3 + i?|e 4 (43) 

with a root 

R L (t) = l/e-eR% (44) 

Accordingly, drops with R(t) > Rl ~ 1/e will grow, lowering the supersaturation e in 
the system. As this proceeds, the monodisperse small droplets (R = Rr{t)) will shrink 
slightly to remain in equilibrium with the current supersaturation. If one neglects this 
last effect, these droplets can play no role, at late times, other than to remove from 
the coarsening process an amount of material uqVb corresponding to that required 
to produce a stable monodisperse trapped emulsion in equilibrium with an infinite 
droplet. The excess material no(Vo — Vb) then becomes concentrated in fewer, larger 
droplets as coarsening proceeds by (essentially) the usual Lifshitz-Slyozov mechanism. 
This scenario is confirmed in the following Section, where we obtain an asymptotic 
solution for the coarsening behaviour. 
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5.2 Asymptotic Analysis 

Following Lifshitz and Slyozov |l| we assume a continuous distribution of droplet sizes 
n(R, t), with n(R, t) dR representing the number density of droplets of radius (R, R + 
dR), at time t. Since droplets of finite size cannot suddenly appear or disappear, the 
conservation of flux through droplet size space requires |TB[ 



dn(R,t) d{n{R,t)U{R,e)) 



Of 



dR 



Conservation of the volume of dispersed phase species requires that 

I f°° Air 

/ R 3 n(R,t)dR + e = e 

Gooff, Jo 6 



(45) 



(46) 



where e = (C — C 00 )/C 00 is the degree of supersaturation and e its initial value. 

Taking Rxit) and -Rl(£) as the two roots of U(R, e) = 0, as defined above, we now 
seek asymptotic solutions of the form 



n(R,t) = A T (R T )f T (R/R T ) + A L {R L )f L {R/R L ) 



(47) 



where fx and fi are two different scaling functions. That is, we assume (see figure |5]) 
that at late times the size distribution splits into two populations, each of fixed shape 
in terms of an appropriate reduced size variable, and each with some time-dependent 
amplitude A. 



n(R,t) 



R 



R 



Figure 5: Approximating n(R,t) by functions centred on R? and R^. 

These two distributions are now assumed to separately obey Eq.^S| for conservation 
of flux through droplet size space. This separation is valid so long as there is no sig- 
nificant range of sizes for which both populations overlap - which is increasingly true 
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at late times. However, the two populations do interact via the supersaturation e. We 
first solve U(Rl,€) = for e(i?i(f)) to obtain 



R L (t) 



R L (ty 



(4£ 



and note that, because of this equation, terms in e and terms in 1/Rl are of the same 
order. In what follows we take the long time limit and therefore write e = 1/Rl + 
0(1/R 3 L ). We then impose conservation of volume fraction, and of the total number of 
droplets (n Q ), to obtain A T and A L . This is done explicitly in Appendix IB] and yields 



An 



n 



R T (e) 



0(e 3 ) 



and 



A, 



e -n {47T/3)R%/(C{oo)v b ) 



(49) 



(50) 



(Air/3)B 0L R L (ty 

where B l is a moment of the scaling function fi defined in Appendix [BJ 

Eq.^5] is then expanded in powers of Rl retaining only the lowest order terms in 
1/Rl (or equivalently, the lowest order terms in powers of e). As shown in Appendix 
P this results in R L (t) = (3^t) 1/3 and 



4/l + Zl 8AMz^ 



dZ L 



7 



h 



■4r + X 



9Ul{Z l )) 



dZ L 



Z L 



+ 0(e) 



(51) 



where Zl = R/RL{t). These equations are identical to those solved by Bray]T8[, re- 
sulting in the consistency requirement 7 = 4/27, and the solution 



h{Z L 



AZ\ cxp 



(3+Z i ) 7 /3(3/2-Z L ) 11 /3 





<Z L < 1.5 
1.5 < Z L 



(52) 



where A is a constant which is determined by the condition / f^Z^dZi = 1. 

The above argument shows that the equation of motion for Rl is exactly the same 
as in a standard coarsening problem (with no trapped species) fl8| , to the leading order 
in small e. Accordingly the solution for Rl{^)i which involves seeking a specific 7 for 
which the scaling distribution remains self-consistent in the long-time limit, is also the 
same. The only difference is in the amplitude Al which, as shown in Eq.^0], has, to the 
leading order in e been shifted by a constant amount corresponding to a reduction in 
eo- In other words, in the long time limit the large droplets behave precisely^Z^ as they 
would for an emulsion with no trapped species but a reduced initial volume fraction 
^0(^0 — Vb)- (Obviously, the latter is assumed positive; otherwise the emulsion is stable 
and will not coarsen.) 



6 Possible Applications 
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6.1 Optimal Sizing of Emulsions 



The use of "trapped" , or less soluble species to stabilize emulsions is widespread in in- 
dustry [P3[ . Such emulsions can be prepared mechanically with various average droplet 
sizes. Let us assume that any trapped species is dissolved at uniform concentration c t 
through the dispersed phase material. As explained in Section |4.4| , to form a stable 
emulsion one must ensure that the initial droplet size is sufficiently large. Roughly 
speaking, this ensures that the typical Laplace pressure 2a /R is smaller than the os- 
motic pressure CtksT. (The quantitative version of this, applicable to any initial size 
distribution, is Eq.|35].) Clearly, it is important to resist the temptation to make the 
initial emulsion too fine (which in ordinary emulsions might be expected to delay coars- 
ening for the maximum possible time). In practice, the need to avoid sedimentation 
(and perhaps coalescence, which we have neglected) may set an optimal initial size close 
to, but above, the absolute stability threshold, Eq.|35]. For the reasons discussed in Sec- 
tion [|, the weaker "spinodal" condition]^] Eq.|2T|, even for a nominally monodisperse 
initial state, cannot guarantee stability. 



6.2 Reversing the Coarsening Process 

Consider a situation where an emulsion is prepared in an unstable state (for example 
by mechanical agitation) and then starts to coarsen. We now ask what will happen 
if, after some time interval, a large number of small droplets (containing a trapped 
species) are added to the emulsion, causing the stability condition V < (V (r])) r) to 
become satisfied for the system as a whole. According to the arguments of Section 



475| , the emulsion is now unconditionally stable. Therefore it will not coarsen further. 
What is more, for a given total volume fraction and a given population of trapped 
species, the final state of the system is the unique one in which droplets of all r] have a 
common chemical potential for the (mobile) disperse-phase species. Thus the addition 
of the small droplets will not only prevent further coarsening, but will in general cause 
previously coarsened droplets to redissolve. Indeed, if there are no trapped species in 
the initial unstable droplets, these will evaporate completely, and their material will be 
entirely absorbed by the added droplets ^] . A formal proof of these remarks, following 



the lines of Section |4.5| , is left to the reader. 

In principle, the coarsening process can be reversed even when it is complete. For 
example, a system of small oil-in-water emulsion droplets with trapped species present 
can, by the evaporation-condensation mechanism, take up oil from an excess bulk 
phase of pure oil. If the stability condition V < (V^r/))^ is met by the system as a 



whole [p6[1 , the bulk phase of oil will disappear entirely (though obviously this process 
may be slow in practice). Like the rest of our conclusions, this one applies only if both 
coalescence, and diffusion of the trapped species through the continuous phase, are 
strictly negligible. With these assumptions, we may minimize the free energy with the 
constraint of fixed trapped species in each drop; the above result for the equilibrium 
state follows immediately. 



6.3 Formation of Mini-emulsions by Shrinking 



"Mini-emulsions" . 



comprised of droplets with radii of between 50 and 150 nm, have 

However their 



many potential uses in industrial and pharmaceutical applications p4 
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formation by traditional mechanical methods, where droplets are formed by strongly 
shearing the ingredients, is limited by the high energy required by the process fl27J, 
and the difficulty of getting a uniform droplet size. An alternative route is to create 
an emulsion of relatively large drops (whose size is also more controllable) and then 
"shrink" them to the required, smaller size. The following describes a non-mechanical 
method which would allow this to be done. 

Consider a situation in which a bulk reservoir of the dispersed-phase species (say, 
oil), containing trapped species at concentration q,, is placed in contact with a stable 
emulsion of oil droplets (in water, say) each containing r\ molecules of the trapped 
species (see figure |]). This system is then allowed to reach equilibrium under the 
evaporation-condensation mechanism. (To speed the process, some gentle agitation of 
the emulsion might be desirable.) 



Trapped and disperse phase 
only. 



Mixture of droplets and 
continuous phase. 



Figure 6: Droplets shrink to attain a stable size, which is in equilibrium with an excess 
bulk phase that also contains a trapped species. 

A final state of equilibrium will be reached when the chemical potential of the 
dispersed phase species in the emulsion equates to that in the bulk. This requires 

nSS = - n L (53) 

where U b fJ^ ,drop are the osmotic pressures of trapped species in the bulk and droplets 
and Ul is the Laplace pressure in the droplets. For ideal solutions this condition reads 

77 2a 
Cb ~ (47r/3)i? 3 ~~ Rk B T 

which is an equation for the final droplet size R. By increasing q, the final droplet size 
can be made as small as one wishes: effectively the oil can be "squeezed out" of the 
emulsion by the osmotic action of the trapped species, now present in the bulk oil phase 
at higher concentration. (This contrasts with the last example in Section |6.2| , where a 
bulk phase, containing no trapped species, could be entirely absorbed by the emulsion 
droplets.) Once the desired size of mini-emulsion is reached, it can be removed and will 
remain stable. 

Note that, because the trapped species is insoluble in the continuous phase (wa- 
ter in this example), there is no need for a semipermeable membrane to prevent its 

21 



(54) 



transfer between the bulk oil phase and the emulsion droplets. Therefore, in practice, 
a more rapid exchange equilibrium might be reached if the "bulk" oil phase instead 
took the form of large macroemulsion droplets which can later be separated out easily 
by sedimentation. Moreover, there is no need for the trapped species in the bulk oil 
phase, and the trapped species in the mini-emulsion droplets, to be identical; so long 
as both are insoluble in water, the same condition for osmotic equilibrium will apply. 
Therefore one can "shrink" an emulsion containing an expensive trapped species (such 
as a fragrance or drug) by contacting it with a cheap polymer solution. 

The shrinking process should not only achieve small droplet sizes, but may also 
allow one to reach concentrations of the trapped phase that would be unattainable by 
normal means. For example, by making an emulsion of dilute polymer solution and 
then shrinking the droplets, it may be possible to achieve within each droplet a highly 
concentrated polymer solution, too viscous to be dispersed mechanically in its own 
right. 

These ideas may be relevant to various encapsulation technologies. Obviously, the 
designation of "oil" and "water" in the above is arbitrary and these could be any two 
phases. We have assumed throughout that coalescence is negligible, which is commonly 
the case for oil-in- water emulsions so long as they contain a surfactant (typically ionic), 
to give a surface repulsion between droplets. The presence of the surfactant should not 
alter our arguments, so long as it is soluble enough in the continuous phase that the 
surface tension o does not vary between droplets. In principle, for small enough mini- 
emulsions, the surfactant could also give rise to significant bending energy terms in the 
free energy of a droplet, which could be included if required (see Section fO] ). 

7 Conclusions 

It has long been known that emulsions containing a sufficiently high concentration 
of a trapped phase will not coarsen via the Lifshitz-SlyozovJTJ mechanism (Ostwald 
Ripening). For practical purposes it is clearly important to know what is "sufficiently 
high" in this context; this information can allow stable emulsions to be designed, rather 
than formulated by trial and error. 

In this work we have provided a general condition (Eq.^) which (in the absence 
of droplet coalescence) will guarantee stability against coarsening, for emulsions of 
arbitrary polydispersity in both initial droplet size and composition. The remarkable 
simplicity of our result stems from a correspondence between the thermodynamics 
of emulsions at fixed droplet number and the equilibrium among multiphase fluids 
(Section [O]). 

Although presented for an ideal solution of trapped species, our condition may be 
generalised to an arbitrary equation of state (see Section [Of ). In fact though, for typ- 
ical parameters (cr ~ 10 _1 N m _1 , R ~ 1/im) stability can be achieved with trapped 
concentrations of order 10 -2 M so that departures from ideality need not arise. An anal- 
ogous condition (Eq. ^0|) was also found that ensures a relatively long-lived emulsion 
even when the "trapped" phase is slightly soluble. 

In deriving Eq.|35| (and its simplified form, Eq.|22] for the monodisperse case) we 
took care to distinguish "nucleated" from "spinodal" coarsening. Eqn.|35] rigorously 
identifies the boundary between fully stable emulsions of arbitrary initial polydispersity 
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(regime I) and emulsions that can lower their free energy by coarsening at fixed droplet 
number. The latter can be subdivided into two classes (regimes II and III). In regime II, 
coarsening can occur by nucleation. The nucleation can arise either because of variations 
in the local environments of droplets, or from the presence in the initial state of even 
a single pre-existing large droplet. We believe that the latter, in particular, can rarely 
be ruled out. 

In contrast, spinodal coarsening (regime III) requires neither mechanism since typ- 
ical droplets are locally unstable. We showed that for monodisperse emulsions, Eq.fHl 
as proposed by Kabalnov et al.JT] actually identifies the onset of spinodal coarsening. 
If in general nucleated processes cannot be ruled out, condition Eq.^l] is insufficient for 
stability. 

The effect of trapped species on the dynamics of unstable, coarsening emulsions 
was also considered. It was found that coarsening would proceed precisely like that 
in an emulsion with no trapped species]]]], but with a reduced initial volume fraction 
n (V - V B ). 

We have emphasised that the emulsions evolution is driven by the competition be- 
tween the osmotic pressure of the trapped species and the Laplace pressure of droplets. 
These principles are exploited in Section |6.3| , where we describe non-mechanical meth- 
ods which might allow the restabilization of a partially coarsened emulsion and the 
formation of stable emulsions consisting of droplets smaller than those attainable by 
traditional mechanical methods. 

There is scope for further work in several directions. For example, it would also be in- 
teresting to consider further the detailed kinetics of emulsions in the metastable regime 
(II); these can be complex, especially if the "trapped" species is itself slightly soluble 
(Section |4.8|) . For practical purposes in designing emulsions, however, this regime is 
best avoided as we have emphasized above. Of greater technological importance is the 
case of emulsions which consist of multiple components of varying solubility. One would 
like criteria for the average droplet composition that will maximise such an emulsion's 
lifetime. We hope to return to this issue in future studies. 
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A Slow coarsening condition for a slightly soluble 
"trapped" species 

In this Appendix we establish Eq.[|0], which is the condition that the size distribution 
remains single-peaked at all times, in a system where the trapped particles are slightly 
soluble. 

We do this by considering the contrary case of a distribution which consists of both 
a rapidly coarsening and a quasi-stable part, with the rapid coarsening of the larger 
drops occuring at a rate determined by the transport of the most soluble component. 
As t — > oo, the size Vl of the larger droplets becomes large, the supersaturation of the 
disperse phases will tend to zero, and the smaller drops will tend to sizes Vb(t)). Here 
r) may now vary among droplets and with time. 
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The total number of droplets is no longer conserved but becomes time-dependent, 
and Eq.|53| is replaced by 

n(0)7(0) - n(t)VE(t) = J n L (r,,t)(VE(t) - V B (r,))d V (55) 

where ni^i], t) is the number of droplets at the larger size Vi(t), and (Vb (£))»? has been 
rewritten as V B (t) for convenience. So if we can ensure that 

m - =gW < (56) 

then 77,^(77, £) = and the distribution must be single-peaked. 
Conservation of the less soluble components requires that 

n(0)rj(0) = n(t)r}(t) (57) 

or 

n(t)V B (t) = v(0)VB(t) m 

n(0) 7?(t) 1 ' 

We next note that Vb(0) = and assume that Vsij]) is convex, as is the case for 
ideal phases where Vb ~ r/ 3 ^ 2 . (This ensures that if 1/2 > 2/1, then 2/1^5(2/2) > 2/2^3(2/1) •) 
If so _ 

T^(t) > Vb(^(t)) (59) 

And 

77(0)V^(t) > mVB(v(t)) > mV B (v(0)) (60) 

To obtain the second inequality, we have again used the fact that V B (rj) increases 
faster than linearly with 77 and also exploited the fact that fj(t) increases with time. 
(The latter follows from the fact that the number of droplets present must decrease as 
time proceeds.) Hence fj(0)V B (t) /r}(i) > Vb(tJ(0)) and 

F(o)-!^ = F(o)-5(||W<F(o)-^ W o)) (6i) 

So if 7(0) < V B {r}{0)) then 7(0) - n(*)T^(*)/n(0) < and n L (r],t) = 0. Hence the 
distribution may only contain a single peak, and will coarsen at rate determined by 
the less soluble component. 



B Determination of At and Al 

In this Appendix, we obtain results for the amplitudes At and Al as defined in Section 



5~2| . These are determined by the constraints of conservation of volume fraction and 
conservation of the total number of droplets. We write n(R,t) = + Ul, with tit = 
A T f T (R/ R T ) and n L = A L f L (R/R L ). We then define fj> and fi to be normalised 
so that J °° fxiZ^dZx = fiiZ^dZi = 1, where Zt = R/Rt and Zl = R/Rl- 
Conservation of total number of droplets then implies 

A T R T + A L R L = n (62) 
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Defining J °° f T Zj,dZ T jC OQ v h = B 0T and / °° fi J Z\dZ L jC (yo v\ ) = B 0L , then conservation 
of volume fraction implies 

47T 

— (A T R 4 T (e)B 0T + A L R 4 L (t)B 0L ) + e = e (63) 

or 

a o M fa - e) - B m A T (^IS)R T (e) t 



Substitution into Eq.|62| for A T then gives 



__ (e - 6) - £? 0T A T (47r/3)iMe) 4 ^ 
T i? T (e) (47T/3) J B Li?L(t) 3 ^T(e) 1 J 

which, to the order required, is Eq.|49|. Likewise for A L we obtain 

A (6 - e) - n B 0T (A7r/3)R T (ef 

Al ~ WJW^aW (66) 

Substituting Eq.g^ for R T (e) and Eq.|4| for e{R L ) yields 

_ 6 - 720^(4^/3)^1 (1 + (3R B /2)B 0T (4n/3)R B ) 
L (4n/3)B 0L R L (t)± (4n/3)B 0L R L (t)s 1 ' 

Since Eq.[H] determines that droplets in the neighbourhood of Rb will relax at an 
exponential rate towards Rt, fr (which is normalized to unity) will converge at an 
exponential rate onto fx = 5(Z T — 1). Hence B 0T equals l/(C(oo)?;&) and A L becomes 

_ 6 - (n (4n/3)R B )/(C(oo)v b ) (1 + (3R B /2)(Air/3)R%)/(C(oo)v b ) 
L (Arc/3)B 0L R L (ty (4tt /3)B QL R L (tf 1 ] 

which to the required order is Eq.[50|. 

C Determination of ul 

We have assumed that riL(R,i) satisfies (by itself) the continuity equation in droplet 
size space, 

dn L (R,t) = d(n L (R,t)U(R,e)) 
dt dR 
Writing n L (R,t) = A L f L (Z L ) we obtain 



(69) 



d(A L f L (Z L )) R Lr7 d(f L (Z L )) dA L 

dt " ~T l ZlAl oz L + -W h{ZL) (70) 

Defining 

K = e - n (4ir/3)R B /C(oo)v b 
1 (47r/3)5 0i 1 ) 

and 

1 + (3R B /2)B 0T (4n/3)R B /C(^)v b 

K2 = wmbJl (72) 
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then we find 



So we obtain 



dA L a Rl a u Rl 



dn L (R,t) R L ( d{f L {Z L )) \ R L 

= ~ Al r- l [ 4h + Zl ^z^~ ) ~ K2 w L h 

We next write d(riL(R,t)U(R, e)) /dR in terms of Z L to obtain 



(73) 



(74) 



R L R 2 L (4f L + Z L 



d(f L (Z L )) - 

dZ L ■ 



h 



dZ L \ Z L 



4{R a /R L ) 2 

(R a / Rl) 2 



Expanding 1/(R? L A L ) to lowest order in 1/Rl we obtain K 2 Rl/R\Al 
keeping only the leading terms on the right hand side we are left with 



R L R 2 L (if L + Z L ^^) + K 2 R L R L l 



d(MZ L )) ( J_ 
Z L 



dZ L 



1 
~7 



(75) 
K 2 R L R L ; 

(76) 



Since the R.H.S. is a constant, then as t — > oo the L.H.S. must also approach a constant 
value. Trying R L = Tt@ with Y a constant, then 



and 



r\r\ = prh 3 ?- 1 



r l r l = prH 2 ^ 1 



(77) 



(78 



Clearly both will tend to zero if (3 < 1/3, which would not correspond to a coarsening 
state. However if f3 > 1/3 then the L.H.S of the equation would diverge at late times. 
Hence (3 = 1/3, and as t — » oo, RlRl — » 0. Writing 7 = T 3 /3 then results in the same 
equation as that obtained and solved by Bray [[18]] , namely Eq.[5l]. 
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